Elastic consequences of a single plastic event : 
a step towards the microscopic modeling of the flow of yield stress fluids. 
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With the eventual aim of describing flowing elasto-plastic materials, we focus on the elementary 
brick of such a flow, a plastic event, and compute the long-range perturbation it elastically induces 
in a medium submitted to a global shear strain. We characterize the effect of a nearby wall on this 
perturbation, and quantify the importance of finite size effects. Although for the sake of simplicity 
most of our explicit formulae deal with a 2D situation, our statements hold for 3D situations as well. 
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I. INTRODUCTION 

An increasing body of experiments on macroscopic flow 
of various complex systems evidence spatially heteroge- 
neous behaviour. We focus here on the (large) sub-class 
of such systems that display a macroscopic yield stress 
(among which foams, suspensions, emulsions, colloidal 
glasses,...). Typically these systems flow homogeneously 
at large stress/shear rate, whereas at intermediate shear 
rate they may exhibit spatial coexistence between a flow- 
ing and a frozen region Q, |^ H, ^ |^ or intermittent het- 
erogeneous flow To this point, there is little insight as 
to whether the mechanisms leading to such macroscopic 
behaviours are generic or dependent on the specific mi- 
croscopic structure of the fluid and the corresponding 
interactions. 

A general class of "elasto-plastic" models has been 
put forward to apprehend those macroscopic behaviour, 
which was first applied to seismologic modelisation [^. 
In these models, the medium first responds elastically to 
a global forcing (either stress or strain). The deforma- 
tion or stress can then locally induce a rearrangement 
or plastic event, if a local threshold is reached. Such a 
plastic event locally relaxes a stress that is elastically re- 
distributed in the medium, and can trigger other local 
events. In this picture, the macroscopic flow is the out- 
come of the collectively organized sequence of local rear- 
rangements. Although this mesoscopic description seems 
very reasonable, many questions remain to be answered 
for this scenario to be operational. First, what is (are) the 
basic plastic event (s), and how can it (they) be identified 
in a given flowing complex material ? Second, what is 
the constitutive (dynamic) equation that describes such 
a single plastic event under a local forcing ? Third, how 
does such an "event", locally relaxing stress, perturb the 
surrounding medium ? The answer to this last question 
is obviously linked to the nature of the plastic event. 

As to the first two questions, i.e. the nature and the 
description of the plastic event, various convincing pic- 
tures have been proposed in the literature. In a pioneer- 



ing work, Bulatov and Argon introduced a phenomeno- 
logical description of a single plastic event, which allowed 
them to describe many properties of macroscopic plastic 
flows |3,B@1- In their simulation, the unit cell can un- 
dergo several fixed plastic deformations specific to their 
hexagonal geometry. Later, on the basis of molecular 
simulations of a Lennard Jones glass under imposed shear 
stress 10], and building on earlier works by Spaepen and 
Argon, Falk and Langer introduced the notion of shear 
transformation zone (STZ), which described a local lim- 
ited zone where rearrangements occur. The occurrence 
of very localized plastic events is most easily evidenced 
in foams [3, where they take the form of Tl rearrange- 
ments. Langer JTI| then constructed an analytical "mean 
field" elasto-plastic model, introducing STZ as zones with 
a plastic tensorial deformation. More recently, Baret et. 
al. T^l, and Braun performed numerical simula- 
tions on lattices, in which a plastic event consists in a 
local scalar displacement occurring when the local stress 
reaches a yield stress. Thus, within a very general class of 
elasto-plastic model, the notion and the description of a 
plastic event is now well documented and clarified. How- 
ever, a clear description of the consequences of a localized 
plastic event on the stress distribution in the material 
(third question) still needs to be constructed. 

This is the purpose of the present paper : using a 
rather general description for a localized plastic event, we 
compute the long-range elastic perturbation that such an 
event induces in an elastic material. We characterize its 
symmetry and amplitude, as well as the way it is mod- 
ified if the event occurs close to a solid boundary. An 
a priori counter-intuitive result which emerges from our 
calculations, is the crucial role played by finite size effects 
in the modeling of flowing elasto-plastic materials. We 
limit ourselves here to the study of the elastic effects of 
a single event, and leave for a later report the analysis of 
the collective organization of the plastic events when the 
material is flowing. 

The paper is organized as follows. In section II we 
specify the general "elasto-plastic" model that we use: we 
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assume that the medium is homogeneous and isotropic, 
as weh as incompressible for simphcity. In section UTTl we 
consider an infinite geometry, and describe a local plastic 
event induced by shearing and the full characterization 
of the perturbation it elastically induces. In this case 
there is no difference between a forcing at imposed stress 
or imposed strain. In section Hvl we focus on finite size 
geometries, where the system is bounded by solid walls. 
First, we describe how a wall attenuates the perturbation 
induced by an event occurring in its vicinity. Secondly, 
we explicit how in a finite size medium, the perturbation 
depends on the global forcing. We calculate the average 
stress relaxation induced by an event at imposed strain, 
and give explicit formulae to compute the corresponding 
stress field relaxation everywhere. In section we con- 
clude and briefly highlight important consequences for 
the modeling of flowing systems. 



sum of a purely elastic response to the forcing (denoted 
with superscripts 0) and of the perturbation induced by 
the occurrence of the plastic event (denoted with super- 
scripts 1). 

With these notations : 



cr = cr° + crl (5) 

where e^"'^ describes the localized plastic event. To go 
any further we need to pay attention to the boundary 
conditions imposed on those fields, which brings us to 
specify the global geometry of the system. We start be- 
low with an infinite medium, before addressing finite-size 
effects in section IV. 



II. ELASTO-PLASTIC MODEL 

We assume, following many of the previously quoted 
studies, that the displacements and deformations are 
given by the simple superimposition of a plastic flow (the 
localized plastic events) and an elastic distortion of the 
medium. We further assume that the medium is homo- 
geneous, isotropic, and linearly elastic. In addition, we 
focus for simplicity on the incompressible case (the com- 
pressible case can be studied following the same lines), 
so that the elastic properties of the medium are fully de- 
scribed by the shear modulus /i. 

Denoting u(r) the total displacement vector at position 
r, the strain tensor is given by e = i(Vu-|- (Vu)*). From 
our hypotheses, this total strain is the sum of an elastic 
strain and a plastic strain (non-zero only at the locus of 
plastic events): 

e^eP^+e"^ (1) 

Incompressibility corresponds to : 

V.M = (2) 

With the hypotheses of linear elasticity and incom- 
pressibility, the total stress tensor is s = —pi + cr where 
p is the pressure and a verifies : 

cr = = 2^e - 2^eP' (3) 

As we have in mind the slow flow of pasty materials 
we neglect inertial effects so that mechanical equilibrium 
simply requires : 

V.(cr-pl) = (4) 

We consider the classical situation where an applied shear 
(either imposed deformation or stress) induces elastic 
loading of the material, up to the point where it triggers 
a single localized plastic event. The consequent state of 
the medium in response to the applied forcing is here the 



III. INFINITE MEDIUM 

We start with the limit case of an infinite medium, 
where a driving at infinity imposes either an applied shear 
strain or an applied shear stress. The purely elastic re- 
sponse e^^'^ is homogeneous. Let us then focus on the 
perturbation 1 generated by a plastic event described by 
e^' . If the system is stress driven, the boundary condition 
for the induced perturbation is: 

cri(oo)^0 (6) 

whereas for a strain controlled system, the boundary con- 
dition reads : 

ui(cx))^0 (7) 

Equations Q , lO and Q can then be reformulated in 
a well defined problem for the perturbation field: 

V.u^ = 

V.(2^ei -/I) = 2^V.eP' (8) 

which must be solved either with cither (jSJ or ((T)). 

Obviously the problem at hand is directly related to 
the response of a purely elastic (incompressible, isotropic, 
homogeneous) system to a punctual force. Let a force f 
act on the medium at the position r'. The displacement 
field u at position r created by this force is given by: 

V.u = 

V.(2/ie -pi) -|-f(5(r) = (9) 

The solution of this system, that actually satisfies simul- 
taneously both types of boundary condition ((t(oo) — > 
and u(oo) — > 0), is the Oseen tensor which is most eas- 
ily dealt with in reciprocal space (we use hats to denote 
Fourier transforms): 

u(q) = 6(q).f 
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FIG. 1: The medium is sheared in the y direction at constant 
strain 7. the infinite system of section III corresponds to 
H 00, U ^ 00 and jj = 7/2 



with: 



6(q) 



(10) 



With this tool we can construct the solution to the 
system This solution is clearly also independent of 
the specific boundary condition (6) or (7), and can be 
simply written: 



ul(q) = 2M6(q).(iq.eP') 



(11) 



At this point, no assumption has been made on the 
nature of the plastic event (fully described here by the 
localized eP'(r), or equivalently by eP'(q)). Therefore 
this expression gives in a very general form the dis- 
placement field induced by a plastic strain in an infinite 
medium (for both type of boundary conditions). By 
derivation and using (3) one easily obtains a similar 
formula in reciprocal space for the stress perturbation 
a^. Both relations lead in real space to expressions for 
the propagators for displacement and stress respectively. 
Rather than producing these formulae in a formal 
general context we specify to a given geometry first. 



Plastic events with a simple shear symmetry in 
2D 



To pursue analytically without dealing with opaque 
tensorial formulae, we now make the simplifying assump- 
tion that the local plastic event has the symmetry of the 
global forcing which we chose to be that of simple shear. 
We further focus on the two dimensional case : the 
H ~^ 00 limit of Fig. ^ We stress however that most 
our conclusions are also valid for the 3D situation (see 
comments further). 

With the hypothesis above, the plastic deformation 

and 



tensor corresponds to simple shear e 



pi 

xy 



^x'xT^yy neglected. Expression ((TT|l becomes: 



pi 



uliq) = 2fj.{d^x-'iqy + Oxy.iqx)^^ 
■Wy(q) = '2'KOxyiqy + Oyy.iq,^)eP^y 



(12) 



For a localized plastic event e^y — eoa^(5(r) (i.e. of 
typical amplitude ep and microscopic spatial extent a^), 
the perturbation displacement described by H12|) is anal- 
ogous to the displacement induced by a set of two dipoles 
of forces (F, 2a) as represented on figureEl with F — afie^ 
(or more precisely its limit for a — > with aF — ^t^a? 
kept constant). 

The shear stress perturbation corresponding to H12|l 
can be obtained using (3), and reads : 



2t^ePjy 



(13) 



Somewhat similar formulae can be obtained for axx and 
ayy but we will in the following mostly focus on the shear 
stress. 

We define formally the propagators P°°, G°° that de- 
scribe the consequences in terms of displacement and 
elastic shear stress of a single plastic event in an infinite 
medium by: 

uHr) = jdv'V^{v-v')ety{v') (14) 

aly{v) = 2^, jdv'G°^{v-v')e%{v') (15) 

Practically, the propagator for the displacement field can 
be derived from equation Eq. (|12|l either with the Fourier 
Transform or by derivation in real space of the Oseen ten- 
sor. The propagator for the shear stress is then deduced 
from linear elasticity. 



FIG. 2: The perturbation due to a localized plastic shear is 
equivalent to the perturbation due to a set of two dipoles of 
forces with F = o/xeo. 
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FIG. 3; Perturbation of the shear stress field for a plastic 
event occurring at the origin in a medium submitted to a 
shear strain or stress. 



In the present two-dimensional geometry, explicit for- 
mulae in reciprocal and real space are: 



G°°(r,6l) 



9 2 



(16) 
(17) 



Hence, in a system forced with a symmetry of simple 
shear, the perturbation of the shear stress due to a local- 
ized plastic event is of quadrupolar symmetry (see Fig. 
and decreases with a power law in two dimensions (in 
three dimensional systems, the quadrupolar symmetry is 
conserved and the power law is -^). 



B. Global effect of a plastic event : 

We now consider more global effects of the same local- 
ized plastic event e^y(r') = a^eoi5(r'), namely quantities 
integrated over a whole layer (of constant y). The fol- 
lowing integral rules are easily derived from the previous 
calculations: 



+ 00 



oo 

+ 00 



-oo 

+ 00 



C7ly{x,y)dx 
ul{x,y)dx 
ul{x,y)dx 



= 



Sign{y)a'^eo 




(18) 



The first equation states that the shear stress resulting 
from the plastic event is redistributed in such a way that 
the integrated stress on every layer is unchanged, i.e. 
there is no net release of stress over a layer, and con- 
sequently no change in the net force applied from above 
on the system! The second relation indicates that the av- 
erage horizontal displacement in a (horizontal) layer de- 
pends only on whether it is above or below the event (but 



not on its distance to the event), while the third equa- 
tion expresses that the average vertical displacement over 
a layer is zero. For simple shear in three dimensions, sim- 
ilar equations hold for quantities integrated over planes 
perpendicular to the loading displacement gradient. 



C. Relation to other studies 

Let us now compare the results we have obtained at the 
end of section A, for the displacement and stress fields 
induced by a single localized plastic event of simple shear 
symmetry, to related studies in the literature. 

First, our results can be compared with the full analyt- 
ical description of an elasto-plastic inclusion in an elastic 
matrix by Eshelby 14*1 . In that study, plastic shear strain 
occurs only within the finite-sized inclusion which yields 
a perturbation of the shear stress around the inclusion. 
Whatever the explicit shape of the inclusion, the long- 
range behaviour of that perturbation is strictly identical 
to the equivalent for the 3D case of expression (|14|) fas we 
have checked). 

To study collective effects, Baret et al. J2] simu- 
lated elements with local yield stresses on a 2D-lattice 
(semi-periodic boundary conditions). In contrast with 
the present study, they modeled the plastic event by a 
simple scalar displacement. As a check, we computed the 
plastic strain tensor corresponding to such a local event 
in expression ((mi , and calculated the corresponding per- 
turbed shear stress. This yields a perturbation of dipolar 
symmetry, consistent with the propagator that they nu- 
merically evaluated on their lattice. This provides a val- 
idation of our procedure, but mostly underlines that the 
nature and symmetry of the elementary plastic event se- 
riously affect the propagator describing its consequences, 
and therefore potentially the collective interplay of such 
events and the resulting macroscopic flow behaviour. We 
believe that the form of plastic event used in the present 
study is more suited for the actual description of the flow 
of elastoplastic materials. 

Kabla and Debregeas (15| performed an explicit nu- 
merical simulation of a two dimensional foam under shear 
strain. In their quasi-static procedure, at each step the 
length of the film is minimized at constant bubble vol- 
ume. They mimic 'Tl' event by reorganizing sets of four 
bubbles when film length decreases below a critical value. 
They studied the stress rearrangements following such 
Tl events in their simulation. Averaging over many such 
events, they found that statistically these stress perturba- 
tions have a quadrupolar symmetry ( with a slight tilt of 
the axes with respect to those of the macroscopic shear 
{x,y), probably due to a structuration of the foam by 
the fiow into a slightly non- isotropic medium) . They ob- 
served that this stress field coincides with that generated 
in an elastic medium by a set of dipoles with the orien- 
tation of the global forcing. What is remarkable is that 
the outcome of their cellular simulation is very consistent 
with the outcome of the (elastic) continuum approach fol- 
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FIG. 4: A shear strain ^ is applied to a two dimensional 
elasto-plastic medium. A plastic event occurs at position r' 
and is here represented by a square, and we seek for the elastic 
deformation of the rest of the medium that will add up on the 
elastic loading symbolized by the parallelograms. 



lowed here. 

To describe the deformation of plastic amorphous ma- 
terials, Langer has introduced the concept of introduces 
shear transformation zones (STZ). In lul, he studies the 
response of a 2D material to an applied deviatoric stress. 
The plastic strain tensor is described without any as- 
sumption on its orientation. The results in that paper 
suggest that for a plastic strain tensor with the symme- 
try of the forcing, the deviatoric stress induced by the 
STZ is of quadrupolar symmetry, with a power law de- 
crease \. Within our model, we find for such an event 
also a quadrupolar symmetry but a power law decrease 
of and do not understand this discrepancy. 

We now turn to the most important part of our paper, 
namely the effect of the finite-size geometry on the stress 
distortion due to a localized event. 



IV. FINITE MEDIUM 

As in the analysis above, we consider here localized 
plastic events with the symmetry of the forcing (simple 
shear), equivalent to a set of two dipoles of forces. There- 
fore the propagator for such a plastic event can in princi- 
ple be viewed as the sum of the propagators for the four 
forces, which brings us back to considerations pertaining 
to the Green's function for the effect of a single force on 
the elastic finite medium. 

We focus on the case of an imposed shear strain rep- 
resented on figure 21 where two solid walls adhering per- 
fectly to the medium are shifted horizontally. 

Therefore the boundary condition for the total dis- 
placement field is: 



u(a;,±i7/2) = ±U 



(19) 



We proceed with the same decomposition as in the infi- 
nite medium case. The displacement field is the sum of 



FIG. 5: The plastic event occurs at the vicinity of the wall : 
r — r'l 3> d. The quadrupolar symmetry of the plastic event 
is assumed not to be affected by the presence of the wall. 



the homogeneous elastic loading and of the perturbation 
due to the plastic event (again we denote by u^,cr^ the 
displacement and deviatoric stress tensor induced by the 
plastic event). The boundary condition for the perturba- 
tion in the imposed strain regime is: 



(20) 



The total response of the elasto-plastic medium is then 



Wi^(a;,y) 
Uy(x,y) 



u. 



H 
1 

y 

IF 



(21) 



Wall effects 



We show here that the perturbation due to a plastic 
event occurring in the vicinity of a wall is more rapidly 
damped than the perturbation due to an elastic event 
occurring in the bulk. 

To be precise, a plastic event is considered to occur in 
the vicinity of the bottom wall (located at yw = —H/2) if 
it occurs at a position r' much closer to the wall than the 
point r where the perturbation is calculated. This con- 
dition requires \H/2 + y'\ ^ |r — r'|, with the notations 
defined on Figure [3 We also focus here on situations 
where the second (top) wall is too far to play a role, that 
is |r — r'l <C — H/2\, so that our approach in this sub- 
section practically deals with a semi-infinite geometry. 

Again, the effect of a plastic event is equivalent to the 
sum of that of the four forces represented on Fig. 
Hence, we first study the displacement field induced by 
a single force F at position r' in the vicinity of the wall. 

The displacement field due to a punctual force F at 
r' in a semi infinite medium is the sum of the displace- 
ment field due to the punctual force and an image with 
respect to the wall in an infinite medium. This com- 
plex image is such that the displacement field it induces 
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FIG. 6: The displacement field due to a punctual force F 
in a semi infinite plane at a position r' is the sum of the 
displacement field the same force would induce at r' in an in- 
finite medium, and that induced by its image with respect to 
the wall also in an infinite medium. The image consists of a 
punctual force — F, two dipoles of forces of strength 2dFx and 
— 2dFy respectively, and two dipoles of potential. The result- 
ing asymptotic long range behaviour is that of a quadrupole 
of strength 2dFx . 



exactly cancel out on the wall the displacement gener- 
ated by the punctual force. Its structure is depicted in 
Fig. Eland recalled below 16]. The location of the im- 
age is the symmetric of the pole r' with respect to the 
wall. This image is the sum of a punctual force — F, two 
dipoles of forces (2(iex, F^ey) and (2dey, Fj,ey), and two 
dipoles of potential (2dex, Fx^/m) and {2dey , — Fyd / fi) 



where d = 



y 



H/2\ is the distance between the event 



and the wall. The analytical expression for the corre- 
sponding Green function is given in |l6j |. Its long-range 
behaviour, |r — r'| 3> d, is dictated by the dipoles of 
forces in Figure 6, yielding a displacement field that 
scales a 



|r-r'|- 



Returning to the effect of a plastic event, where the 
source is now a set of dipoles as in Figure 2, one could 
expect a similar cancellation of the first order terms, and 
thus a far-field displacement scaling as j^^pjj. However 
inspection of the structure of the image show that such is 
not the case (see Figure[7|l : the dipole of horizontal forces 
is duplicated in its direct image, and the set of dipoles 
of vertical forces it generates yield a dipole of strength 
AaF. The contributions of the (original) dipole of verti- 
cal forces are weaker due to cancelations. Altogether one 
is left with a dominant term that has the same geometry 
than in an infinite medium (and a similar l/r decay for 
the displacement) with an amplitude twice as strong. We 
have analytically checked that the stress decay is consis- 
tent with the above picture, as for a localized event close 
to the wall we obtain a propagator that is directly related 
to its equivalent in the absence of the wall Eq. lfT7|) : 



Gwan(r, 0) = 2G°°(r, 0)(1 -f 0(d/r)) 



(22) 



where clearly x — x' = rcos{9) , y — y' — rsm{9). This 
picture is also consistent with numerical calculations of 
the propagator for an event next to the wall in a finite 
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FIG. 7: The displacement field due to two force dipoles in a 
semi infinite plane is the sum of (i) the displacement field it 
would induce in an infinite medium and (ii) the displacement 
field due to its image with respect to the wall. This image 
can be constructed by summing those of the four forces. The 
contributions represented here (not at scale) are the direct 
images (solid arrows, all of amplitude F) and the induced 
force dipoles (dashed arrows, of amplitude 2{d—a)F and 2{d+ 
a)F for the top and bottom one, of amplitude 2dF for the 
two others). Inspection shows that the net dominant long- 
range behaviour is that of a quadrupole twice as strong as 
the original one 



FIG. 8: Perturbation of the shear stress field for a plastic 
event (plastic deformation of amplitude eo = 1) occurring 
next to the top wall. 



geometry, which yields the picture in Fig. |S1 (calculation 
to be described in subsection IVB below). 



B. Propagator in a finite thickness medium 

In this subsection, we turn to a medium of finite thick- 
ness H . We focus again on an imposed strain situation. 
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FIG. 9: A periodic array of plastic events at position x' + 
kL,y'. 



and therefore on the problem corresponding to the sys- 
tem of equation (jH)) together with the no displacement 
boundary conditions (pn|l . We indicate ways of calculat- 
ing the propagator in this geometry but mostly empha- 
size consequences of a single event on integral quantities. 



1. Finite H , Infinite L 

A first method to treat the case of a medium of finite 
thickness and infinite length, consists in the systematic 
construction of a series of images so as to cancel the dis- 
placements on both walls. 

In the previous subsection we reported the explicit 
structure of the image of a punctual force needed to can- 
cel the displacement due to this force on one of the walls. 
The image with respect to the wall at -ff/2, unfortunately 
creates a displacement at the wall — 7?/2, so that its own 
image with respect to the wall —H/2 has to be consid- 
ered, and this process has to be repeated for every image. 
Thus formally two infinite sums of images are required 
to express the displacement for a punctual force within 
two walls. Following this strategy, Pozrikidis per- 
formed a full analytical calculation of the deformation 
field induced by a point force. 

Formally, the perturbation due to a plastic event can 
then be deduced by adding up the consequences of each 
of the four forces it consists of. This leads to expressions 
that although exact are heavy to deal with and some- 
what opaque. We therefore turn to other methods in 
the following, focusing on a geometry periodic in the x 
direction. 



2. Finite H , L Periodic : first method 

We now focus on a system of thickness H and of fi- 
nite extent L in the x direction, and consider periodic 
boundary conditions in that direction. This is equiva- 
lent to analyzing in an unbounded geometry the effect 
of a periodic array of plastic events of a given ampli- 
tude at positions {x' + kL,y') with fc G Z (Figure (SJ, or 
formally i^y{r) = a'^<^oJ2k&^i^ ~ + kLx)). The re- 
sulting displacement field can then be viewed as the 
sum of the displacement field induced by a similar peri- 
odic array of plastic events in an infinitely thick medium 



u°° and of a correction term due to the finite size H. 
— ^joa _j_ function of x — x' , y and y' . 

The displacement field u°° induced by the periodic ar- 
ray of plastic events in an infinite medium can be ex- 
pressed using the propagator for a single event H14() : 



(23) 



Given its periodicity in the x direction, this expression 
can be formally decomposed in Fourier series: 

oc 

u-(a:, y) = U-(y) + ^ ( U- (y) cosi2mTx/L) 

n=l 

+ lJZ{y)sm{2mTx/L)) 

(24) 

From equation 1)181) . the components of the zcroth-mode 
vector function are: 



Sign {y-y')a^eo/L 

(25) 



The correction displacement field is the solution of 
incompressible linear elasticity with no source (i.e. the 
set of equations H3I4() with no plasticity) but with the 
boundary conditions required to grant that is zero on 
the walls: 

v"ix,±H/2) + u°°{x,±H/2)^0 (26) 
is thus also periodic and can be written : 



v"ix,y)^Y^{y) + Y, i^cniv) cos{2nnx/L) 

Vfjy) sm{2mTx/L)) 



n=l 



where the functions V^, can be calculated inde- 
pendently (i.e mode by mode) with the boundary con- 
ditions: Vf„(±i?/2) = -U-(±iJ/2), Vf„(±if/2) - 
— U^j(±iJ/2). We skip the full calculation in this subsec- 
tion (which can be performed e. g. a s in the low Reynolds 
number hydrodynamic study in ^^), as we display in the 
next subsection an exact expression of the displacement 
calculated in a framework that is more amenable to nu- 
merical simulation. 

Instead we focus here on the zeroth mode V^(y). 
From equations (|25|l and (|26|l . its components are simply: 



-2a^eo — 



0. 



(27) 



This suffices to deduce consequences of the plastic 
event in terms of integrals over constant y lines: 



u'dx = L{V^iy) + lJ^{y)) 

J-L/2 



8 



so that 

J-L/2 ^ 
.L/2 

/ ul{x,y)dx = 

J-L/2 

Then, hnear elasticity implies that the overall variation 
of the force on the top plate due to the plastic event is: 

SF ^ a^y'dx = -2^ (28) 

J-L/2 

The corresponding drop of the average shear stress in the 
medium is obviously S{a) = SF/L. 

This exact expression shows that a single plastic event 
(of a given fixed amplitude t^a^ ) results in the release 
of the net force exerted by the medium on the walls by 
a quantity scaling as jj. Remarkably, this quantity is 
independent on the position of the event in the medium. 
A corollary is that a finite density of plastic events 0, 
will relax this total force by an amount ^ {(j)HL)6F, 
corresponding to a relaxation of the average stress S{a) = 
6F/L independent of the size of the system. 

Note that the integrals over a period L calculated 
above for the consequences of a periodic array of plas- 
tic events, are equal to the integrals over x from —oo to 
+CX3 in a finite H infinite L geometry for the consequences 
of a single plastic event. For example, in the latter ge- 
ometry SF = f^^cTxy^dx — —2 ^'^'^ , which clarifies in 
what sense equation (|18|l obtained in the previous section 
for an infinite medium corresponds to the limit H ^ oo. 
The global relaxation due to a single plastic event is con- 
sequently directly related to the finite size of the system. 

An alternative derivation of this finite-size dependence 
of the force (or average stress) relaxation is proposed in 
appendix which yields a somewhat complementary 
physical insight. 

3. Finite H , L Periodic : second method 

Although we aim to focus in this paper on the quali- 
tative aspects presented previously, we propose here an 
explicit formula for the stress perturbation induced by 
a localized plastic event in a finite size geometry. The 
actual expression is presented in reciprocal space, which 
may appear at first somewhat cumbersome, but turns 
out to be convenient in numerical calculations. For sake 
of readability, we provide here only the principles of 
this two step derivation and the resulting formulae, the 
details of the calculations being presented in appendix IbI 

First step: we formally extend the actual system 
{0 < y < H, and L periodic along x) by the follow- 
ing two operations: first an antisymmetric image sys- 
tem is constructed that extends to —H < y < 0, then 
the 2H thick resulting system is repeated periodically 




-15 -10 -5 5 ID 15 



FIG. 10: Perturbation of the shear stress field for a plastic 
event (plastic deformation of amplitude eoa'^ = 1) occurring 
in the vicinity the wall. The discretization corresponds to 
H/a = L/a = 32 

in the y direction. If the plastic strain in the orig- 
inal system is eP''{x,y), our construction yields a sys- 
tem without walls that is 2H periodic in the y direc- 
tion, with in the upper half period H > y > a plas- 
tic strain eP''*{x,y) — eP'(a;,j/), and in the lower half 
—H < y < an antisymmetric image plastic strain 
e^'*(a;,y) = — eP'(a;, — y). This construction ensures that 
the overall displacement generated by these strains has 
an x component that is symmetric by reflection by the 
planes y = and y = H, and an y component that is 
antisymmetric in the same operations. We have there- 
fore generated a solution that satisfies the condition of a 
zero y component of the displacement on the y — and 
y = H planes (loci of the walls in the original system) . 

Second step: We now want to cancel the remaining dis- 
placements along X without modifying the above result, 
and without adding sources in the system < y < H . 
This can be achieved by adding on the planes y — 
and y — H appropriate force fields fx directed along x 
(again asymmetric and 2H periodic along y). Given the 
symmetry and periodicity of the system it is clear that 
the y displacement on the walls is not modified by this 
addition. 

When this is achieved, we have in the upper half < 
y < H, a, solution to JSJ that satisfies the no displacement 
boundary condition on the walls H20|l . The corresponding 
stress field can be expressed in Fourier series : 

a\x,y)= J2 e'f"^e*«"^(T*(m,n) 

with p„i ~ and <?„ — . It is the sum of the term 
directly generated by the plastic strains e* and that due 
to the added force fields: 

a*{m,n) ^ 2^l{G°°im,n)iP^*{m,n) 

(m, n) + iqnOxx{m, n))fx{m, n)} 

o-*(0,0) = (29) 
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The propagator G°° and the Oseen tensor O are the di- 
rect counterparts for periodic systems of those defined 
in section Hill Of course the force field fx in the above 
expression is itself proportional to the plastic strain; the 
corresponding formulae are given in appendix IbI 

Inverting back to real space allows to compute numer- 
ically the response to a localized event. We have checked 
that for a large system the stress created by a plastic 
event far from the walls according to this formula is close 
to that obtained analytically in section III for an infinite 
medium. Similarly, for an event directly neighbour to 
the wall, we recover the results of subsection IVA for a 
semi-infinite medium. Fig. llOl represents an intermediate 
situation : the event occurs in the vicinity of the wall. 

V. CONCLUSIONS AND PERSPECTIVES 

Starting from a general elasto-plastic model, we have 
computed in different 2D geometries the modification of 
the shear stress resulting from a localized plastic event 
with a symmetry of simple shear. We have first calcu- 
lated the corresponding perturbation in an infinite sys- 
tem forced with a symmetry of shear. The stress field 
is of quadrupolar symmetry and decreases with a power 
law 4j. in two dimensions. Then, we showed that the 
stress field perturbation due to a plastic event occurring 
close to a wall has a modified near field structure but de- 
cays far away with the same law and pattern, although 
with an amplitude twice as large. Eventually, we have 
proposed two ways of calculating the perturbation of the 
stress field due to a plastic event occurring in a finite 
medium. The first one allowed us to demonstrate in a 
simple way that a plastic event of a given strain ampli- 
tude relaxes the average stress by an amount which is 
independent of its position (i.e. distance to the walls) 
and inversely proportional to the size of the system. The 
second one allowed us to derive explicit expressions that 
permits calculation of the whole stress field in a finite-size 
geometry. 

The extension of our results to a three dimensional 
situation is rather straightforward, and the qualitative 



statements are obviously similar. Also we have focused 
on a scalar description of plastic strain and induced 
stress, but the same steps can be taken if one seeks to de- 
scribe the whole tensorial stress field generated by plas- 
tic events of arbitrary symmetry. Eventually, we have 
focused in section IV on a situation where the walls were 
kept fixed during the plastic event. The extension of our 
results to situations where the force on the plates is kept 
fixed is immediate: one simply needs to add to our so- 
lution a simple shear displacement (corresponding to a 
stress —5F/L with the notations of section IV). 

This study is meant to be a first step towards the 
modeling of the flow of an elasto-plastic material. The 
next one consists in plugging in a plastic law that de- 
scribes the onset and evolution of the localized plastic 
events. Coupling such a local plastic behavior to the 
long-range elasticity described here should yield interest- 
ing collective behaviours and hopefully insights in the 
flow mechanisms. For such an endeavour our quantifica- 
tion of finite-size effects is important for the steady-state 
average balance between the stress released by the plas- 
tic events and that imposed by the elastic loading. Also 
the geometry and decay law of the elastic perturbation 
is to be considered when addressing the emergence of 
a collective/cooperative organization at low shear rates. 
Important questions regarding possible spatial and tem- 
poral heterogeneities in such flows (0, 0) S IE E| could 
be addressed at the light of the present results: 

- We have shown that a localized plastic event relaxes 
the average stress but also modifies the stress pattern on 
all the lines of the system decreasing the stress of some 
elements but increasing that of others. Does this allow 
shear banding with a limited zone fiowing in coexistence 
with a non flowing region ? 

- The average stress released during an event of given 
strain amplitude is independent of the position of the 
event in the medium, and is the same on all lines. Yet 
the geometry of the propagator is clearly modified by the 
proximity of a wall. Do these elements favor a localiza- 
tion of the fiow at the wall ? 

These questions will be addressed in a forthcoming pub- 
lication. 
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APPENDIX A: FINITE SIZE EFFECTS. 

We propose an alternative simple approach which gives 
a different physical insight into the effects of finite size ev- 
idenced subsection IIV Bl The following argument holds 
equally for an infinite L — > oo geometry or for an L pe- 
riodic medium, and is here presented for the latter. We 
consider a plastic event occurring at (cc', y') described by 
(^y{x, y) — eoa^S{x — x' ,y — y') (and its repeated images 
along X due to the L periodicity). The stress perturba- 
tion is: 

.L/2 

'^lyix^y)= dxi dyiG"^{x - xi,y,yi)ePy{xi,yi) 

J-L/2 J-H/2 

(Al) 

which actually defines the propagator for the 

present periodic medium. The corresponding force re- 
lease on the bottom of a layer at height y is 

rL/2 .L/2 

6F= dx aly{x,y) ^ aho dx G""^ {x - x' ,y,y') 

J-L/2 J-L/2 

(A2) 

We remark that this can be rewritten 

nL/2 nH/2 

SF= / dxi / dyiG"^{x - xi,y,yi){eoa'^6{yi -y')) 

J-L/2 J-H/2 

(A3) 

which means (see (Al)) that the variation of the inte- 
grated stress over a layer at height y {5F) induced by an 
event occurring at position a;',?/', is equal to the stress 
variation at a site (x', y) induced by a continuous sum of 
events on a layer y' . 

Representing these events by pairs of dipoles, the sum 
reduces to that of the x components as the y components 
cancel out upon summation on the line, so that we are 
looking for the stress generated at {x', y) by the two lines 
of forces aX y' ~ a and y' + a depicted on Figure ITTI 
The displacement field resulting from each of these sums 
separately are easily calculated and depicted on the right 
of Figure ^] For example the second one is the solution 
of the following problem: on each side of the plane y = 
y' — a, the medium is purely elastic, with no displacement 
at the walls, and at y = y' — a a, continuous displacement 
and a jump in the stress of amplitude /leoa. 

The total effect of the continuous line of events is the 
displacement field presented on the left of Figure^] The 
slopes of the profile for y > y' + a, and y < y' — a are 

2 

the same, and correspond to SF — —2 ^^ ^° , as derived 
differently in the main text, see equation H28|) . Remark- 
ably, this value, which corresponds to the force release on 
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FIG. 11: Left: a continuous line of events at height y' is 
equivalent, using the representation of Figure 2 for a single 
event, to two lines of horizontal forces a,t y' + a and y' — 
a. Right: displacement profiles resulting from the sum of 
negative forces at y' — a (left), and from the sum of positive 
forces at height y' + a (right). 



the top wall, is independent of the position of the plastic 
event. 

u u 




FIG. 12: Left : displacement profile due to the line of events 
at y' (obtained by summing the two profiles of Fig. Illl right'). 
and schematic description of the elastic loading. Right: total 
response to the global forcing and a line of plastic event (a 
homogeneous fracture), obtained by summing the two profiles 
on the left. 

We remark that the deformation profile obtained 
here is very similar to the one obtained by Kabla and 
Debregeas in their explicit simulation of a shear foam 
(which in their case was the line average consequences of 
a single local event, equivalent to the present local con- 
sequence of a line of events according to the argument 
below (A3)). 

As a follow-up along the same picture, the cumulative 
displacement after loading plus a continuous line of plas- 
tic events is represented on the right of Figure El : the 
total shear in the system is then the homogeneous shear 
due to the forcing minus the shear release due to the 
plastic events. 

APPENDIX B: STRESS FIELD INDUCED BY A 
PLASTIC DEFORMATION IN A FINITE 
MEDIUM 

In this appendix we provide the explicit derivation of 
the stress field change due to a plastic activity in a finite 
system, follow the strategy outlined in subsection IVB.3. 
The aim is to obtain the response of the system depicted 
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on figure 1^1 to a strain solicitation e^'(x, j/) with "stick" 
boundary conditions on the wall similar to|2Dl For this 
calculation we choose to take the origin of the y axis on 
the bottom wall so that the walls (and thus the boundary 
conditions) are at y = and y = H. 



From our calculations for an infinite medium in section 
mil it is logical to write the total displacement using 
propagators for plastic events and the Oseen propagator 
for simple forces: 
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FIG. 13: Auxiliary system corresponding to the real system 
of figure IHl it is periodic in both direction and symmetric in 
the y direction. The driving plastic strain field e^'* is anti- 
symmetric. 

a. Geometry of the auxiliary system- To solve this 
problem, an auxiliary system with no walls is con- 
structed. Its geometry is that of the original system, 
symmetrized with respect to y = 0, and periodized along 
y. The new system has thus a period L in the x di- 
rection, and a period 2H in the y direction, see Fig. 
US Fields in this new system are described by stars. 
The symmetry with respect to the plane y = is made 
in such a way that the plastic strain is antisymmetric 
eP'*(a;, —y) = —eP''*{x, y), which is equivalent in the force 
dipole representations to having the force fields in the 
lower half to be the symmetric of that in the upper part 
F*(x,y) = F*{x, —y). The symmetry and periodicity of 
the problem ensures that the displacement field u* gen- 
erated by this plastic strain satisfies verify by symmetry 
and periodicity : 



u*Jx,Q)^Q-ul{x,H)=Q 



(Bl) 



h. Response to a plastic event- This new elastic 
system is submitted to a plastic strain eP^*{x,y) = 
eP'{x,y) if y > 0, and eP'*(a;,y) = -eP'(a;, -y) if y < 0. 
For reasons explained in the text (section IVB.3), a force 
field /* is added at the location of the walls in the real 
system, to cancel the displacements Ux{x, 0) and Uxix, H) 
created by the plastic strains. 

/:(x, y) = r^{x)5{y) + rH{x)5{y - H) (B2) 

This "no displacement on planes y = and y = W con- 
dition is implemented in reciprocal space. We introduce 
notations for the Fourier series: 



M*(x,y) 

u*(m, n) 
2'Km 



1 



2HL 



dxdy u* (x, y) e^' ~ 



27rn 



q 



P-n 



u*{m,n) = P°^{m,n)iP'-*{m,n) 

+ d{m,n)Um,n) (B4) 



The equations for P°° and the Oseen tensor O in the 
present representation (Fourier series) are formally ex- 
actly similar to those obtained in @ and H14|) in terms 
of Fourier transforms (with g^, and qy — > </„). 

The condition u*{x,0) = 0:u*{x,H) = allows to 
compute f*{m,n). 



fxim,n) 



for odd n 

E2 fp;,(g.^-p;^)g^'*(m,nO 

for even n and m ^ 
for even n 



(B3) 



^•j^ describes the sum over odd values of n', and ^2 f^e 
sum over even values of n' . 

The shear stress in this geometry can then be derived 
from the displacement field: 



&*{m,n) = 2fi{^{iprnul{rn,n) + iqnU*{rri,n)-eP'-*{rn,n)) 

(B5) 

For symmetry reasons, cr*(0, 0) — 0. 

c. Back to the real system - The displacement field 
u* (x, y) that we have constructed and computed satisfies 
the elasticity equations for > y > with the plastic 
strain eP'(x,y) as only source in that area, and verifies 
the boundary conditions u*{x,0) — and u*{x,H) — 0. 
It is therefore the solution to our initial problem for the 
real system: for iJ > y > 0, u-'-(x, y) = u*(x, y), and 
a^{x,y) = a*{x,y). Hence, the expression of the stress 
as given in Equation (|29|l . 



